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INTRODUCTION 

In a quantum simulation experiment, the behavior of a com¬ 
plex quantum model is examined using controllable system, 
which acts as the simulator [1-3]. Collections of trapped 
atomic ions have emerged as excellent standards for quan¬ 
tum simulations of interacting spin models. These electrically 
charged particles are confined in an electromagnetic trap, mi¬ 
crometers apart from one another [4, 5]. Thanks to precise 
spin manipulation capabilities, near-perfect spin readout, and 
a variety of cooling [6, 7] and dynamical decoupling tech¬ 
niques [8-16], trapped ions can be made to follow target 
Hamiltonians with high fidelity [17, 18], making them one of 
the most promising candidates for quantum simulation. 

To date, trapped ions have mostly been used to simulate 
spin one-half Hamiltonians, showing the phase transition from 
the (anti)ferromagnetic to paramagnetic phases in the Ising 
model [19-30] and long range correlation functions in the XX 
model [31, 32]. By enlarging the spin’s degree and moving 
into integer spin chains, new and subtle physics can appear 
[33-41]; for example, the local orders vanish and we are left 
with hidden orders only [42]. 

In last few decades considerable efforts have been made 
to investigate the non-local physics which appears in integer 
spin chains, e.g., a spin-one XXZ antiferromagnetic (AFM) 
Hamiltonian 

H = ^7,., + +7)^ • (D 

i-j i 

In 1983, Haldane conjectured [35] that Heisenberg chains 
of integer spins with nearest neighbors antiferromagnetic in¬ 
teractions are gapped, unlike the gap-less half-integer spin 
chains. This energy gap corresponds to short-range exponen¬ 
tially decaying correlation functions, compared to the long- 
range power law decaying correlations of half-integer spin 
systems. Later, Den Nijs and Rommelse [36] suggested that 
the Haldane phase of the spin one chain is governed by a hid¬ 
den order, which can be characterized by a nonlocal string 
order parameter. It is consistent with a full breaking of a hid¬ 
den Z 2 X Z 2 symmetry, which was revealed using a nonlocal 
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Figure 1. The phase diagram of the spin-one XXZ antiferromag¬ 
netic (AFM) Hamiltonian (eq. 1). The first derivation in this text 
(as well as the previous derivation [40]) cover 0 < A < 2 (in blue), 
whereas the second derivation in this text covers the whole positive 
A > 0 plane (in blue and red). 


unitary transformation by Kennedy and Tasaki [37, 38]. In 
2010 Pollmann et al. [39] showed that the Haldane phase 
can also be described by the doubly degenerate entanglement 
spectrum. These characteristics hint that the Haldane phase 
is a topologically protected phase in one dimension. Hence, 
quantum simulation of the Haldane phase enables the explo¬ 
ration of topological behavior in relatively simple systems. 

Recently an experiment by Senko et al. [41] has used 
trapped ions to simulate spins in higher degrees. In this ex¬ 
periment a spin-one Hamiltonian was engineered using state- 
dependent laser forces, paving the way toward highly control¬ 
lable quantum simulation experiments that exhibit hidden or¬ 
ders. Previously, we proposed how to engineer the Hamilto¬ 
nian of Eq. 1 using microwave-based forces on the trapped 
ion spins [40]. In this previous derivation of the Hamiltonian, 
we could cover only 0 < A < 2 of the phase diagram (Fig. 
1). A natural way to extrapolate the microwave approach to a 
laser based implementation is to use non co-propagating Ra¬ 
man beams instead of a magnetic field gradient [43] to induce 
the spin-phonon coupling. 
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In this paper, we first propose a new approach to achieving 
the same result as in [40], but with a laser based implemen¬ 
tation. Then, we propose a new approach for quantum engi¬ 
neering eq. 1, while covering the whole positive A > 0 plane 
of the phase diagram (Fig. 1). Instead of having the XX (flip- 
flop) Hamiltonian SySy -I- S‘^Si as our starting point, in the new 
approach we engineer a Mplmer-Sprensen (MS) like gate [44] 
for spin-one systems that takes the form of an Ising Hamilto¬ 
nian s[si. 

We have found three ways to implement the MS gate for 
spin-one systems [44^6]. These schemes were originally 
proposed for two-level system (qubits), but can be extrapo¬ 
lated to three-level systems (qutrits). For simplicity, we con¬ 
sider the first and second approaches for laser based designs, 
although a similar derivation exists for the third one, and also 
for microwave based implementations. 

We explain how to move experimentally to the interaction 
picture; namely, how to measure in the appropriate basis. Fur¬ 
thermore, we give a detailed explanation of the adiabatic path 
for reaching the Haldane phase. We stress that our model is 
robust to the main noise sources in experiments, exploiting 
dressing fields as a dynamical decoupling technique, and op¬ 
erating in a decoherence free subspace, enabling longer adia¬ 
batic evolution times. 


MODEL 


The phonons — In our model we have N trapped ions each 
of mass ni and electric charge e, forming a linear chain along 
the z axis, The ion formation is determined according to 
the equilibrium of the Coulomb repulsion between the ions, 
and the trapping forces. The vibrations of the ions around their 
settled points, Ar,H, are solved in the harmonic approximation 
obtaining the normal modes and v“, which are the eigen¬ 
states and the eigenvalues of the mode and the ion in the 
a direction respectively. Thus, the ion displacements are rep¬ 
resented as Aria = Ln ^Jhjlmv^ {bla+bna^, and the vi¬ 
bration Hamiltonian is represented as Hyn, = Y.n.a 
where b]ia, bna are the creation and annihilation of a phonon 
respectively [47-49]. 

The spin — The second quantum degree of freedom in this 
model is the spin. For our derivation, we consider the hyper- 
fine strucmre of the ^^*Yb+ion (fig. 2), with microwave en¬ 
ergy separation between the singlet and the three triplet states 
[50, 51]. Removing the mf = 0 state (|0')) from the triplet, 
we are left with three energy levels for modeling the simulated 
spin-one system. Quantum manipulation can occur using Ra¬ 
man transitions via virtual excited states, or using microwave 
driving fields. Note that any other ion having three different 
energy states in the microwave regime would suffice for mod¬ 
eling the spin-one particle, as long as there is a virtual excited 
state, through which the Raman transitions between the three 
levels can occur [52]. 

The interaction — In trapped ion systems, the spin-spin in- 
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Figure 2. The *^*Yb+ground states used for modeling the spin- 
one particle (in green). The beat frequencies between non-co- 
propagating Raman beams (in red) drive the 5 detuned transitions 
between | —1) —>■ |0) and |1) —> |0) (in blue), generating the spin-laser 
interaction in eq. 2). 


teraction takes place by exchanging a virtual phonon between 
the different ions; therefore, a spin-phonon coupling term is 
needed. Using laser Raman beams, the spin-phonon coupling 
term is achieved by a sufficiently large Lamb-Dicke parame¬ 
ter of the non co-propagating beams [52]. For antiferromag¬ 
netic interactions, we typically choose a radial vibration mode 
as the mediator between the spins [53, 54]. Therefore, three 
counter-propagating Raman beams having a momentum dif¬ 
ference along the radial direction perform the two beat fre¬ 
quencies corresponding to the 5 detuned transitions between 
| —1) —> |0) and |1) |0), with a 7t phase difference (fig. 2). 

In the interaction picture with respect to the radial vibration 
Hamiltonian and the bare state energy structure, we obtain the 
red sideband transition for a three-level system; 



After applying the Lamb-Dicke approximation 
Vjj^ns/Nn +1 <C 1, where Nn is the average number of 
phonons of the mode, r\jjx = yjh /2ot v„ is the 

Lamb-Dicke parameter, ki is the laser’s wavenumber, and 
fI_ = s/2 ^|1)-' (0| -I- lO)-' (—1|^ . In the second-order pertur¬ 
bation approach, this results in a XX Hamiltonian in addition 
to a residual term [41]; 


Hxx = Y^jf\ {FiFj +f;fJ) (1 - dij) - ^^5ij I , 
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Our model can be decoupled from the residual term if 
is approximately uniform across the chain, such that Hres 
Y,iFz- model’s relevant ground 

states, as we will see later in this paper, have a vanishing 
projection over the z axis, thus the experiment is performed 
in the decoherence-free subspace. As was suggested in Ref. 
[41], this residual term can be eliminated by adding additional 
beat frequencies that generate the blue sideband transitions to¬ 
gether with the red sideband ones, such that the MS transitions 
for spin-one systems is realized. 

Moreover, using trapped ions, the spin-spin interaction is 
not the nearest-neighbor-only as the classic Haldane work has 
considered, but rather a power-law decay [53, 54]. However, 
as Ref [56, 57] shows, the Haldane phase can still be found 
in power-law interacting systems. 

Up to this stage of the derivation of the AFM XXZ Hamil¬ 
tonian, we have shown how to generate the first two terms of 
eq. 1 . We now pursue the two other tunable terms, the D and 
X terms. 

Generation of the control anisotropy D-term — Generat¬ 
ing the tunable D-term can be done in two ways: (1) using 
an additional transition to generate the A.C. Stark shift, from 
which the D-term will arise (Fig. 3 a); (2) imposing detun¬ 
ings on all the previous driving fields,such that the D-term is 
set aside from the bare state energy structure (Fig. 3 b). To 
generate the A.C. Stark shift we can apply an additional Ao 
detuned microwave driving field, corresponding to the transi¬ 
tion between the two zero states of the triplet and the singlet, 
|0') ^—> |0) (fig. 3 a). Similarly we can obtain the same result 
by applying co-propagating Raman beams corresponding to 
the same transition. Thus we obtain: 
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Figure 3. Generating the £>-term Generating the D-term in eq. 1 is 
done in two ways: (1) (a.) Applying two additional co-propagating 
Raman beams, corresponding to a detuned transition between 

|0') t—>■ |0) results in an A.C. Stark shift ^ |0) (0|. This could al¬ 
ternatively be generated using a Ao detuned microwave driving field. 
(2) (b.) Within the non-copropagating Raman beams (red) (fig. 2), 
we regard the |0) state as D'-shifted, and move to the interaction pic¬ 
ture with respect to the D'-shifted bare energy structure. Thus, we are 
left with —D' |0) (0|, which is O'F^. In our later derivation (which 
covers the entire A > 0 plane), we will impose a 2D' detuning rather 
than a D' detuning. 
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^\0){0'\e'^^' +h.c, (5) 

which yields an A.C. Stark shift term D'fFff- with D' = — 
in the second order perturbation approach, assuming 
Ad- 

Generating a tunable D-term can alternatively be achieved 
by regarding the |0) bare state energy level as D'-shifted, cor¬ 
responding to its real value in the hyperfine structure (fig. 3 
b). In other words, we impose D'-detuning on all the driv¬ 
ing fields generating the XX Hamiltonian (eq. 3). Therefore, 
we are left with the term —D' |0) (0| = D'{Fz)^. In order to 
tune this term we have to shift all the driving fields transitions 
continuously, as was recently shown in [41]. Note that the 
anisotropy D-term is slightly different from these D' terms, 
since the model is 9 rotated while generating the A-term, as is 
explained below. 


Figure 4. Generating D'F)cos6 for the A -generating trick. There 
are two ways of generating this term: (1) (a.) Additional copropagat¬ 
ing Raman beams (red) perform the iAj detuned transitions between 
|±1) -fl |0) respectively, with a Rabi frequency of £1~. These effec¬ 
tive transitions can alternatively be realized using detuned microwave 
driving fields (blue). Therefore, we obtain the desired term from the 
A.C. Stark shifts where D'cosO = G^/Aj. (2) (b.) Adding detun¬ 
ings to the non-copropagating Raman beams (red) (fig. 2), such that 
5' = 5 + Q.'cos 0. In this way we are left with the desired term, af¬ 
ter moving to the interaction picture with respect to the detuned bare 
energy structure. 

Generation of the Ising-like X-term — The Ising-like 
anisotropy term XJijF/F/ is produced using a technique, 
which is revealed by considering the interaction picture. It 
is done by adding a spin operator term, 0-rotated from the z 
axis, namely, QlF^p = Q! {F^cos 9 + Fa sin 9), where a can 
be either or y or their superposition (XY plane), and gener- 
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Figures. Generating fl'fa sin 6 for the A -generating trick, (a.) 
We apply carrier transitions using additional copropagating Raman 
beams (red), or alternatively microwave driving fields (blue), while 
keeping the imposed detunings for generating the D-term (fig. 3 a) 
and the Q.'F^cos 6 term (fig. 4 a). Flere, Hear = H'sin6y/2, and 
a = Cti — 0!_i is determined by the initial phase difference of the 
two effective transitions. The same carrier transitions, only without 
detunings, are used in our later approach as dressing fields that trans¬ 
form the system into the dressed state basis (b.), thus protecting it 
from the magnetic noise. 


We have obtained the required spin-one XXZ AFM Hamilto¬ 
nian, where we can only cover 0 < A < 2. Experimentally, 
the following parameters are realistic; V„ = 2n ■ 5MHz, 5 = 
2n ■ 5.1MHz, Q.^2n- 500kHz,T 7 „j « 0.14, Q! = 2n- lOkHz, 
which gives ~ 2% ■ IkHz, such that Q <D <2%- lOkHz. 

An important advantage over the previous derivation in 
[40], is that in order to obtain the antiferromagnetic interac¬ 
tions, 5 should be larger than the secular frequencies of the ra¬ 
dial vibrational modes, rather than the Rabi frequency n/v/2. 
Since one of the main experimental challenges is to reach a 
high Rabi frequency, in the present derivation we overcome 
this technical obstacle. As for the residual terms we dropped 
on the way, such as the carrier transition that yields a four- 
photon A.C. Stark shift and was neglected from eq. 2, and 
the neglected fast rotating terms and the residual from eq. 3, 
these terms operate with Since the model’s relevant 

ground states have a vanishing projection over the z axis, we 
are decoupled from these undesired terms. 

In this scheme, we move more than once to different inter¬ 
action pictures. In the following section, we explain how this 
can be done experimentally. 


ating each component separately. cos 9 can be produced 
similarly to the two ways that the D term was generated. The 
first way is to drive the transitions between |±1) F-)- |0) with 
±A^ detunings respectively. These transitions can be gener¬ 
ated using co-propagating Raman beams, or alternatively us¬ 
ing microwave driving fields, with a Rabi frequency such 
that the effective A.C. Stark shifts result in the Q/F^cosO 
term, with O! cos 9 = (Fig. 4 a). The second way to 

induce H'/f. cos 9 is by imposing detunings on the non co¬ 
propagating Raman beams; i.e., by imposing ±n' cos 9 detun¬ 
ings on the transitions between |±1) -H> |0) respectively, such 
that the cos 9 term is set aside from the bare state energy 
structure (Fig. 4 b). 

The second term n'Fasin© is obtained using the carrier 
transitions, where their relative phase difference determines 
the resulting operator direction a. The carrier transitions can 
be applied either using co propagating laser Raman beams, or 
simply using microwave driving fields, with a Rabi frequency 
Hear — n'sin0v/2 (fig. 5 a). 

By applying a 9 rotation around the perpendicular axis, 
such that F^ 0 is transformed into F^, the new spin operators 
and the new basis are 9 rotated as well. If we move to the 
interaction picture with respect to the new term we have built 

H'F^, and use the RWA where we require ^ <C «C 

5 — V„, we end up with the following effective Hamiltonian: 
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MOVING TO THE INTERACTION PICTURE 
EXPERIMENTALLY 

Since the measurement is taken in the lab frame and the 
effective Hamiltonian is derived after moving to several inter¬ 
action pictures, we now show how the dynamics of the state 
in the lab frame relates to the dynamics determined by the last 
interaction picture. For this purpose, we take the general case 
where in the lab frame the evolution of the Schrodinger state 
is described as follows; 


Us{t) = exp {-iHs{t')dt') (7) 

|'P,(f))=t/.(f)l%(0)) 

where T is the time ordering operator. Moving to the interac¬ 
tion picture, the dynamics is described as; 

U^{t) = exp (///of) 

Hj^{t) = U^{t)Hi„,{t)Uo{t) 

ft ( 8 ) 

Uioit) = T J^exp{-iHjg{t')dt') 

|'P,„(f))=t//„(f)|'P/„(0)), 

Usually, the first interaction picture Hq is time-independent 
since it is the bare state energy structure. However, in 
general, Hq can also be time-dependent, where //^(f) = 
exp (/ jQHo{d)ddy Here, we only assume that Hq is a single 
particle operator or a global rotation; i.e., it does not create en¬ 
tanglement. In fact, if the realization of the D control parame¬ 
ter is achieved using detunings, varying it during the Haldane 
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phase simulation experiment effectively results in having the 
first interaction picture //o(f) time-dependent. 

Using the definition of the interaction picture 

|‘P/o(0)=t/o(f)l‘P.(0), (9) 

we obtain the relation between the interaction and the 
Schrodinger pictures: 

U,,(t)=ul{t)Us{t), (10) 

since |'P/o(0)) = |'P,(0)). 

Suppose that the first interaction frame is by itself a 
’’Schrodinger” frame for the next interaction picture, such 
that: 

( 11 ) 

where H\ (t) is a single particle operator in the first interaction 
frame. As before 

uj (f) = exp J Hi {t')dt'^ 

t//j {t) = T J exp {—iHi^ (f')) dt' 

thus, 

Ui^{t) = Ui\t)Ui^[t). (13) 

Substituting eq. 10 in eq. 13 we obtain 

t/,(f) = t/o(f)t/i(f)t//,(f); |'P,(f)) =t/o(f)t/i(f)|'P/,(f)) 

(14) 

For A + 1 interaction pictures in our derivation, we obtain 

U,{t) = UQ{t)Ui{t)...UN{t)Ui,{t) 

|%(f))=t/o(f)t/i(f)-t/iv(OI'I'A(0)- 

In order to move experimentally to the interaction picture, 
such that the Schrodinger state will evolve according to the 
last interaction Hamiltonian, we have to rotate the system back 
to counter Uo{t)Ui{t)...UN{t). There are two main ways to 
move experimentally to the interaction picture, as will be dis¬ 
cussed next. 

The first and straightforward way is by applying N concate¬ 
nated global rotations for each interaction picture, except for 
the first one. Suppose we want to measure the system after 
time T, in which all the driving fields were on. In order to 
counter the last interaction unitary Un{t), we only block the 
driving fields that generate the simulated Hamiltonian Ui^{'c), 
while operating with all the driving fields that are responsible 
for the interaction pictures. This will be done for time Tat, such 
that we obtain: 

Us{i + 1 ^) = Uo{t+ + (16) 


and tM is determined by + = I. The same approach 

can be used to counter all the other interaction unitary, obtain¬ 
ing the following pulse sequence: 

(1) . All the driving fields for the experiment time T. 

(2) . All the driving fields that are responsible for the Hi (t) to 
HN{t) interaction pictures for fjv time. 

(3) . All the driving fields that are responsible for the Hi (t) to 
HN-iit) interaction pictures forfAi_i time. 

(N). Only the driving fields that are responsible for the Hi (f) 
interaction pictures for ti time. 

In that way, we obtain the following relation between the 
last interaction picture and the Schrodinger one: 

Us{t + ti\l + ...+ti) = 

Uo{t + + ... + ti)Ul{T + In + ... + ti)...UN{T + t7v)t//^ (t). 

(17) 

Setting t4(T -f fiv + ... + 4) = I for all 1 < k < A, the 
Schrodinger state therefore evolves according to the last in¬ 
teraction picture 

Us{t + tN +... +ti) = Uo{t + tN +... +ti)Ui^{T:). (18) 

with an additional phase Uq{t ++ ... -l-fi), which we do 
not measure. If we want to measure in any different basis, we 
can rotate the system using any one of 14 (t -f -f ... -f 4) = 

The second way to move experimentally to the interaction 
picture is simpler. Since the pre-factor of 14 ^ (t) is a global ro¬ 
tation t/o(T)t/i (t)... 14 v(t) = we can apply just one 

global rotation to counter the total interaction unitary rotations 
all together. Since any rotation in the Bloch sphere can be gen¬ 
erated by two independent rotations (e.g. two orthogonal rota¬ 
tions), it is sufficient to use the first and second interaction pic¬ 
tures, namely Ho{t) and Hi{t), which are orthogonal rotations 
for that task. Measuring in any other basis can be achieved 
by an additional rotation, which can be added to all the inter¬ 
action unitary, and thus be represented effectively with H(,{t) 
and Hi{t). Next, we explicitly show how to apply the first 
approach in the two derivations of the XXZ-D Hamiltonian. 

Experimental realization — While generating the XXZ-D 
Hamiltonian we move twice to different interaction pictures: 
the first one is with respect to the bare energy gap of the qutrit, 
yielding the following interaction unitary: 

t/o(T) = exp(-/[(a)o-£>')|0) (0|+I1 11) (1| - |-1) (-1|] t) 

(19) 

where we used the detuning approach to generate the D 
term. The second interaction picture is with respect to the 
A—generating term, resulting in the following interaction uni¬ 
tary: 

t/i (t) = exp (—/n' [cos dF^ sin OFa] t) , (20) 

with a = x,y. Moving to these interaction pictures results in 
the effective Hamiltonian (eq. 6) which yields the simulated 
evolution t//j (t). 
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As was discussed above, in order to move to the interaction 
frame, we first have to shut down the effective Hamiltonian. 
This is done by blocking the two non-copropagating Raman 
beams (setting Q. = 0), and eliminating the D' term. Shut¬ 
ting down the effective Hamiltonian should take fi such that 
Ui{x + ti) = 1 . 

There is no need to counter Uq since it only yields an un¬ 
measured phase. Leaving the experiment at this point results 
in measuring the state in the basis. Now, if we want to mea¬ 
sure in the simulated basis, namely = cos 9F^ + sin OFa for 
a = x,y, we need to rotate the system around a perpendic¬ 
ular axis. Since by using the carrier transitions we generate 
Q! sin OFa, we have to rotate the system with an Fp operation, 
with p ± a. This can be achieved after countering t/i, by a 
% 12 phase change of the laser or microwave driving field, such 
that instead of operating with Fx the driving field will operate 
with Fy. In a similar way we can measure in any basis we 
choose. 


ADIABATIC QUANTUM SIMULATION 

Once we know how to quantum engineer the Hamiltonian 
of the system under investigation with tunable parameters, the 
adiabatic quantum simulation proceeds as explained below. 
We initialize the system in a trivial phase by appropriately set¬ 
ting the parameters of the Hamiltonian, and we initialize the 
system’s state in its trivial ground state. Then, we change the 
Hamiltonian parameters adiabatically, slower than the energy 
gap, such that the system stays in the ground state of the in¬ 
stantaneous Hamiltonian, until it reaches the ground state of 
the non-trivial phase that we want to investigate. 

In our case the system under investigation is the Antiferro¬ 
magnetic XXZ-D spin-1 Hamiltonian, and its trivial phase is 
the large-D phase, with a tensor product of |0) in each site, as 
its trivial ground state. Since we are operating with a rotated 
basis F^ = cos OF^ + sin OFa, we can initialize the system in a 
tensor product of |0) in theT’j. basis using polarization [50, 51]. 
Then we can rotate the state with orthogonal operations: Fp 
(with a _L P), similarly to what was described above. Now, 
we set a large D parameter of the Hamiltonian, and start the 
adiabatic variation of the D parameter until we reach the Hal¬ 
dane phase regime. 

When getting closer to the thermodynamic limit, the en¬ 
ergy gap in the second order phase transition closes [58]; thus 
the adiabatic approximation cannot hold for increasingly long 
chains. Overcoming this problem takes advantage of the fact 
that the Haldane phase is a symmetry-protected topological 
phase. Thus, we can use a symmetry-breaking perturbation 
in order to go around the D ^ H phase transition, while still 
operating adiabatically (Fig. 6). 

The Haldane phase is protected by the following symme¬ 
tries: a bond centered spatial inversion Sj S^j+\, a time 
reversal symmetry Sj —> —Sj and the dihedral D 2 symme¬ 
try, which is a ;r rotation around the x, y and z axes. In or¬ 
der to break all these symmetries we add a perturbation term 



Figure 6. Adiabatic path using the symmetry breaking perturba¬ 
tion. When crossing a second order phase transition the energy gap 
is closed, and the adiabatic approximation is invalid. We break the 
symmetries of the Hamiltonian (eq. 1) and thus, go around the phase 
transition, keeping a finite energy gap during the whole path. 

Hpert = Since we know how to engineer the 

term, which is F^ in our derivation, the symmetry break¬ 
ing term can be produced by individual addressing, which is 
achieved by focusing the laser beams. 

Note that since the quantum simulation experiment is adia¬ 
batic, its duration has a lower bound determined by the energy 
gap. However, like any other quantum experiment, it also has 
an upper bound that is determined by the coherence time. If 
the experiment lasts longer than the coherence time, decoher¬ 
ence processes might destroy the quantum information carried 
by the system, and the results cannot be trusted. Hence, im¬ 
munity to the main noise sources is crucial. 


ROBUSTNESS OF THE GROUND STATES TO NOISE 

To benefit from a long coherence time we need to be de¬ 
coupled from the main noise sources. The most fidelity dam¬ 
aging noise source is the ambient magnetic field. Usually, to 
simulate spin-1/2 systems, the ^^*Yb+clock states are used. 
However, when simulating spin-1, the use of the Zeeman lev¬ 
els leaves the system vulnerable to magnetic field fluctuations. 
The second most important noise source is the fluctuations of 
the Rabi frequencies of the non-copropagating Raman beams 
that generate the spin-spin interaction £1, the carrier transition 
Clear, and the driving fields that generate n'sin0Fa. In or¬ 
der to counter these noise sources, we combine the continuous 
version of the dynamical decoupling technique and the unique 
quality of this model that makes it possible to conduct the ex¬ 
periment in a decoherence-free subspace. In our scheme, we 
use driving fields that refocus the noise in directions perpen¬ 
dicular to the final basis we operate with (i/.). Specifically, 
these driving fields operate as dressing fields performing dy¬ 
namical decoupling, thus, we are left with noise sources in 
the z direction. However, since the relevant ground states of 
the large D and the Haldane phases belong to the decoherence 
free subspace, we are protected against these noise terms. 

The ground state of the large D phase is the topologically 
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trivial state of a tensor product of |0), and the ground state of 
the Haldane phase [37, 38] has the same number of sites occu¬ 
pied by 11) and | — 1) , where 11) | O) and | — 1) are the eigen¬ 
states of with eigenvalues 1,0, — 1 respectively. Therefore, 

these ground states are eigenstates of with a zero eigen¬ 

value, such that our model is decoupled from this operation 
and thus from these noise sources. Thus, the quantum simula¬ 
tion operates in the decoherence free subspace. For the same 
reason, all the neglected A.C. Stark shifts resulting in Y.jFz 
could have been dropped in the derivation. 


VERIFICATION OF THE HALDANE PHASE’S GROUND 
STATE 


The ground states of the Haldane phase are characterized 
by; (1) an excitation gap and exponentially decaying correla¬ 
tions of the local order C^{i — j) = (S'g,S a), where () denotes 


the expectation value in the ground state, (2) a nonvanishing 


where 

is the string corre- 


nonlocal string order (//) = lim 

ii-j) = {-S'a exp in 1 S'a 

lation function, (3) a symmetry-protected double degenerate 
entanglement spectrum, obtained by dividing the systems into 
two parts, tracing out one of them and diagonalizing the re¬ 
duced density matrix [59]. 


Experimentally verifying the ground states of the Haldane 
phase can be accomplished by directly measuring the correla¬ 
tion functions of the local orders and the string order (signa¬ 
tures (1) and (2)). For simplicity, suppose we want to mea¬ 
sure Cj-{i — j) or equivalently — j). Experimentally we 
can only measure the dark singlet state |0) without the ability 
to distinguish between the bright triplet states | ± 1), yet by 
single addressing we can rotate a chosen spin and measure the 
other states as well. Since (s[S{\ = Pj’-j +P'_^\ —P‘{Li — 


p',j 

■^-1.1 


,where is the probability to measure <^|a)'(a| \b)^ 
four measurements would suffice [60]. The same holds for the 
string correlations, except here we have to count the bright 




c/ 

+ 1 


. To 


triplet states of the intermediate spins, exp 
measure the correlations in other basis, we first globally rotate 
the system to the desired basis, and then implement the same 
procedure. 


A NEW APPROACH FOR COVERING THE WHOLE 
POSITIVE A > 0 PLANE 


In the above scheme, the Ising-like X parameter is limited 
0 < A < 2, and we can not cover the whole phase diagram, 
just like in the derivation of the previous paper [40]. This has 
to do with the term, where a _L j3, in the XX Hamil¬ 
tonian (eq.3). In order to solve it, we have to suppress this 
limiting term; namely, we have to perform a MS like gate of 
spin-one, rather than the red side-band transition which results 
in the limiting XX Hamiltonian. Once we generate the 
interaction term instead of the XX Hamiltonian, we can span 
the whole positive A plane of the phase diagram by using the 
same A—generating trick as above. 

As was suggested by Senko et al [41], the straight-forward 
way to generate the MS transitions is by applying additional 
beat frequencies to drive the blue sideband transitions in addi¬ 
tion to the red sideband ones, based on Ref [44]. Yet, in this 
section, we would like to show another way to realize the MS 
Hamiltonian, based on the Bermudez et al. gate proposal [45]. 
Here, we use the carrier transitions that were used to generate 
n'sindAa in the above derivation, in order to suppress the 
limiting term, and generate the F^Fx interaction (Ising Hamil¬ 
tonian). Thus, in the rotating frame of the bare state energy 
levels, the carrier transitions yield 


_ y. j 

^^car — 7 


( 21 ) 


where a is in the XY plane, and is determined by the rela¬ 
tive phase of the two carrier transitions (fig. 5 b). For sim¬ 
plicity, we assume that a = x corresponding to a vanishing 
initial phase between the carrier transitions. These transitions 
dress our qutrits, and have significant implications in terms of 
suppressing the magnetic field noise. Moving to the dressed 
state basis can be thought of as a — ;r/2 rotation about the y 
axis; thus, Fx,Fy,F^ in the bare state basis are transformed to 
S^^,Sy,—Sx in the dressed state basis, respectively. In the ro¬ 
tating frame of the dressed state energy structure (eq. 21), the 
red sideband transition (eq. 2) becomes: 


F^red — 






^‘■car t 
s/l 



^‘■car t 
v/2 



{b]^e"'"* +h.c). 

( 22 ) 


If 5 — ^ Q.carj'JX we can neglect the fast rotating terms; 

thus we are left with the MS Hamiltonian for qutrits: 



Full tomography is simply infeasible experimentally for in¬ 
creasingly long chains. However, in order to measure the en¬ 
tanglement spectrum (signature (3)), we can make tomogra¬ 
phy of a part of the system only. Namely, we can measure the 
reduced density matrix of this part containing a few spins, and 
diagonalizing it numerically. 


resulting in the Ising Hamiltonian for spin-one systems, in the 
second order perturbation theory, if „/2v/2 ^ 5 — V„: 


Hi, 


Sing 


i<j 


/ 

S[s{{l-5,j) + 

V 



(24) 
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This looks like we are not proceeding with the derivation of 
eq. 1, but rather are going backwards. We have already had 
the XX Hamiltonian (eq. 3), and with additional effort (eq. 
21), we only obtain the Ising Hamiltonian. However, using 
the trick for generating the A—like term, we also generate the 
XXZ Hamiltonian, this time covering the whole positive X 
phase diagram, unlike the previous derivation. In the follow¬ 
ing, we pursue a tunable D term, and a tunable A term. 

Generation of the control anisotropy D-term — As before, 
generating the tunable D-term can be done in two ways; gen¬ 
erating an A.C. Stark shift, or imposing detunings. To gener¬ 
ate the detuned transitions to implement the first method we 
can apply an additional Ao detuned microwave driving field, 
corresponding to the transition between the two zero states of 
the triplet and the singlet, |0') <—^ |0). In a similar way, we 
can obtain the same result by applying co-propagating Raman 
beams corresponding to the same transition (fig. 3 a). We then 
follow the previous steps of derivation; by first transforming 
to the dressed state basis, where |0) —^ (|m) — |t/)) /v^, and 
then moving to the rotating frame of we obtain 

^ A) - W (tfl >') +h.c. 

(25) 

In the second perturbation theory, assuming ^ 

D? r~ 

A/), and <C v2Q.car, only the A.C. Stark shifts survive, 
whereas the other off-diagonal terms coming from the Ra¬ 
man transitions between \u) < —> \d) are suppressed in the 
RWA. Thus, we are left with the desired term, setting 

D' — ^ 

Generating a tunable D-term can alternatively be achieved 
by regarding the |0) bare state energy level as 2D'-shifted cor¬ 
responding to its real value in the hyperfine structure (fig. 3 
b). That is to say, we impose 2D'-detuning on all the driv¬ 
ing fields generating the Ising Hamiltonian (eq. 24). There¬ 
fore, we are left with the term 2D' |0) (0|. Transforming to 
the dressed-state basis and moving to the interaction picture 
with respect to as was mentioned above, all the off- 

diagonal terms are suppressed using the RWA, if we assume 
that D' ^ \fl£lcar- Once again, we are left with DfSff'. 

Generation of the Ising-like X-term — To implement 
the A-generating trick; i.e., adding a spin operator term 
0-rotated from the dressed z axis, which is = 

n'(5jCos 0 4-5c[ sin 0), we generate each component sepa¬ 
rately. The first cos 0 term is easily produced since it 
can be set aside from the carrier transition (eq. 21), such 
that instead of eq. 21, — O'cos 0^ 5^ is used as the 

dressing term. Regarding the second term D'5asin0 there 
are two alternatives. Choosing a = y, the term n'5vSin0 
is produced by using additional microwave driving fields, 
corresponding to the two transitions between |=f1) o |0) 
with ± ^ — Q! cos 0^ detunings, ± | initial phases and the 

same Rabi frequency Q.y = s/lQ! sin 0 respectively (fig. 7 a). 
Similarly, we can obtain the same result with co-propagating 



Figure 7. Generating 0!Sa sin 9 for the A-generating trick. There 
are two ways to realize this term of the A-generating trick; (1) 
(a.) to engineer il'Sy sin 6, we apply additional co-propagating laser 
Raman beams, which perform the two Raman transitions between 
1^1) O |0) with ±5;^ = ± — fl'cosdj detunings, initial 

phases and the same effective Rabi frequency fij = \/2D' sin 9 re¬ 
spectively. (2) (b.) to engineer —D'sin05;c we apply a z-polarized 
radio-wave driving field, on resonance with the dressed-state energy 
structure (fig. 5 b), and with a Rabi frequency of £1^ = 2Q! sin 6. Re¬ 
gardless of method, we obtain the desired term by moving to the 
rotating frame of the dressed energy structure, and using RWA while 
assuming O! < VlQ.. 


laser Raman beams corresponding to these transitions. 

Choosing a = x, the —Sx sin 0 term, is simply obtained by 
applying an additional z-polarized radio-frequency (RF) driv¬ 
ing field, cos ^ ^— Q! cos ^> with Q.^ = 2Q! sin 0 
(fig. 7 b). Thus, by moving to the rotating frame of the bare- 
state energy structure, the RF driving field is not affected. 
Finally, for both a = x,y, the desired terms are obtained in 
the rotating frame corresponding to —Q'cosoJ S^, us¬ 
ing the RWA where we assume that Q! <C \/2Q.car- 

By applying a 0 rotation around the perpendicular axis, 
such that 5^ 0 is transformed into S^, the new spin operators 
and the double dressed state basis are 0 rotated as well. If 
we move to the interaction picture with respect to the new 
term we have built Q!S^ and use the RWA where we require 
(n?? • 

<C n' <C 5 — V„, we end up with the following effec¬ 
tive Hamiltonian; 





_ _, 


Sy^x T 


cos^ 0 


-f5'5^sin2 0 


re//' 

il 


I 


/ COS^ 0 


-sin^O^ ( 502 . 


(26) 


Therefore, we have obtained the required spin-one XXZ AFM 
Hamiltonian, while covering the whole positive A > 0 plane in 
the phase diagram. Experimentally, the following parameters 
are realistic; Vn=2n ■ 5MHz, d =2n ■ 5. IMHz, Q.car = 2% ■ 
IMHz, D! = 2n- lOkHz, D. = 27Z- 500kHz ,ri„j « 0.14,which 
gives 7*/^ = 27Z ■ IkHz , such that 0 < D < 2;rl0kHz. 
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EXPEIUMENTAL REALIZATION OF MOVING TO THE 
INTERACTION PICTURE 

In this scheme of generating the XXZ-D Hamiltonian, we 
move to three interaction pictures: the first one is with respect 
to the bare energy gap of the qutrit, which yields the same first 
interaction unitary Uo{t) (eq. 19) of the first approach. The 
second interaction picture is with respect to the microwave 
dressed state energy gap, giving rise to the following interac¬ 
tion unitary: 

(27) 

and the last interaction picture is with respect to a superposi¬ 
tion of the 0—rotated term used for generating the X— gener¬ 
ating trick. It results in the following interaction unitary: 

U 2 {'u) = exp (—/n' [cos OFx -I- sin OFa] t) , (28) 

with a = z,y. Moving to these interaction pictures results in 
the effective Hamiltonian (eq. 26), which yields the simulated 
evolution U[^{x). 

As was discussed above, in order to move to the interaction 
frame, we first have to shut down the effective Hamiltonian. 
This is done by blocking the two non-copropagating Raman 
beams (setting Q. = 0), and eliminating the D' term for t 2 time 
duration, such that 1 / 2 ( 7 + f 2 ) = I- 

The next stage is countering Ui. For that purpose, we have 
to shut down the driving fields responsible for 1 / 2 - Specifically 
we have to block the driving fields that generate Q! sin OFa, 
and to change the Rabi frequency of the carrier transitions 
to Q.car/— Q.’cos 9. This stage should take fi, such that 
Ui{t + t 2 +ti) =1. As before, there is no need to counter Uq, 
and we can measure in any basis we desire. 

Similarly to the above scheme, during the adiabatic quan¬ 
tum simulation experiment, we counter the main noise sources 
- the ambient magnetic field fluctuations, and the Rabi fre¬ 
quency fluctuations - using a combination of the continuous 
dynamical decoupling technique with decoherence-free sub¬ 
space. 


Ui{t) — exp ( —i 


^car 


— Q ! cos 0 


SUMMARY 

We have discussed how to quantum engineer the spin-one 
XXZ-D AFM Hamiltonian with two schemes. The first cov¬ 
ers 0 < A < 2, whereas the second covers the whole positive 
A > 0 plane in the phase diagram. It enables us to explore the 
regime of the Heisenberg Hamiltonian of integer spin systems 
where the Haldane phase resides. We have explained how to 
adiabatically generate this non-trivial topological phase, start¬ 
ing from the large D phase with its trivial ground state. During 
the adiabatic path, the ground states are robust to the fluctua¬ 
tions in the magnetic field and Rabi frequencies, and belong 
to a decoherence free subspace. This permits a longer adia¬ 
batic path with higher fidelities. We have also shown how the 


Haldane phase can be verified with simple experimental mea¬ 
surements. This proposal may thus constitute an important 
step towards exploring topological phases with trapped ions. 
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